Single-site Anderson Model. I Diagrammatic theory 
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The diagrammatic theory is proposed for the strongly correlated impurity Anderson model. The 
strongly correlated impurity electrons are hybridized with free conduction electrons. For this system 
the new diagrammatic approach is formulated. The linked cluster theorem for vacuum diagrams 
is proved and the Dyson type equations for electron propagators of both electron subsystems are 
established, together with such equations for mixed propagators. The approximations based on 
the summing the infinite series of diagrams are proposed, which close the system of equations and 
permit the investigation of the system's properties. 
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I. INTRODUCTION 

The study of strongly-correlated electron systems be- 
come in the last decade one of the most active fields 
of condensed matter physics. The properties of these 
systems can not be described by Fermi liquid theory. 
One of the important models of strongly correlated elec- 
trons is the single-site or impurity model introduced by 
Anderson' 1 ! in the 1961 and discussed intensively in a 
lot of papers' 2-15 ' . It is a model for a system of free 
conduction electrons that interact with the system of lo- 
cal spin, treated as just another electrons of d- or /- 
shells of an impurity atom. The impurity electrons are 
strongly correlated because of strong Coulomb repulsion 
and they undergo the exchange and hybridization with 
conduction electrons. This model has some properties 
similar to those of Kondo model having more interesting 
physics' 16-18 !. It has the application for heavy fermion 
systems where the local impurity orbital is / - orbital. In- 
vestigations of impurity Anderson model have used inten- 
sively the methods and results obtained for Kondo model 
by Nagaoka' 18 ' and other authors' 19 ' 20 ] . All the cited pa- 
pers are based on the method of equation of motions for 
retarded and advanced quantum Green's functions pro- 
posed by Bogoliubov and Tiablikov' 21 ' and developed in 
papers' 22-24 ' . 

The first attempt to develop the diagrammatic theory 
for this problem was realized in the paper' 25 !. These au- 
thors used the expansion by cumulants for averages of 
products of Hubbard transfer operators and their alge- 
bra. 

With introduction of Dynamical Mean Field Theory 
the interest for Anderson impurity model increases be- 
cause infinite dimensional lattice models can be mapped 
onto effective impurity models together with a self- 
consistency condition' 26,27 !. 

The Hamiltonian of the model is written as 

H = Hq + H in t , 



H 



pr c _i_ rrf 



n o 

Hint 
„ f 



5>(k)C+C kCT . 

e / ^/+/ CT + C/nfn{ ! (1) 

(7 



N 

fa fo i 



where Ck^C^) and f a (fa) ~ annihilation (creation) op- 
erators of conduction and impurity electrons with spin a 
correspondingly. e(k) is the kinetic energy of the con- 
duction band state (k, a) , e/ is the local energy of / - 
electrons, U - is the on-site Coulomb repulsion of the im- 
purity electrons and N is the number of lattice sites. Hi nt 
is the hybridization interaction between conduction and 
localized electrons. Summation over k will be changed to 
an integral over the energy e(k) with the density of state 
po(e) of conduction electrons and the matrix elements will 
be considered as the function of energy V(e). Because of 
the hybridization term of the Hamiltonian down and up 
spins of conduction electrons come and go in the local or- 
bital and there is no appearance of spin flip process. Thus 
the important parameters of the Anderson model are the 
band width W, the conduction density of states p(e), the 
local site energy e/ and the on-site Coulomb interaction 
U. The electron energies are counted of chemical poten- 
tial jU of the system: e(k) = £(k) — jU, e/ = e/ — . 
There is also an energy parameter T(e) associated with 
the hybridization term 

r ( £ ) = jf E v ^ ~ £ ( k )) = ^ 2 ( £ Mc)- (2) 



This function is assumed to be a constant, independent of 
energy. The term in the Hamiltonian involving U comes 
from on-site Coulomb interaction between two impurity 
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electrons. U it is far to large to be treated by pertur- 
bation theory. It must be included in Ho which is non 
interacting Hamiltonian. The existence of this term in- 
validates Wick's theorem for local electrons. Therefore, 
first of all, we formulate the generalized Wick's theorem 
(GWT) for local electrons, preserving the ordinary Wick 
theorem for conduction electrons. Our GWT really is the 
identity which determines the irreducible Green's func- 
tions or Kubo cumulants. Such definitions have already 
been used by us for discussing the properties of one-band 
Hubbard model^ 28-30 ! and the formulation of the new 
diagram technique for itl 31 ~ 34 l. 

In Section II, we start by introducing the temperature 
Green's functions for the conduction and impurity elec- 
trons in interaction representation, formulate the gen- 
eralized Wick theorem and provide explicit examples of 
diagram calculation for thermodynamical potential and 
full propagators. The results are analyzed in Section III 
and compared to the other data in Section IV. Some ap- 
proximations are discussed in Section V and in Section 
VI there are the conclusions. 



tion have the form: 



G(k,<7,r|k>,V) = - (TC kff (r)C kv (r')[/(/?)};, 
g(a,r\a',T') = (T/ ff (r)/ CT ,(r')C/(/3))o • (3) 



Besides them there are also anomalous ones: 



F(k,cr, t | -k, — a' , t') = - (T<7 kff (r)<7_ k ,_ CT , (r')C/(/3))o , 
F(-k,-a,r\k',a',T') = - (TC_ k _ CT (r)C k v'(r'M/?))o , 
f(a,r\-a',r') = (T/ ff (r)/_ CT , (r')t/(/?))o , 
7(-<t,t | <j',t') = (T/_ CT (r)7 CT ,(r')C/(/3))o, 

(4) 



II. DIAGRAMMATICAL THEORY 

The Matsubara renormalizcd Green's functions of con- 
duction and impurity electrons in interaction representa- 



if the system is in superconducting state. Here r and r' 
stand for imaginary time with < r < (3, j3 - inverse 
temperature and T is the chronological ordering opera- 
tor. The evolution operator U (/3) is given by 
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U{(3) = Texp(- f H mt { T )dT) = Y, { —^f- j dn- f dT n T(H int (T 1 )...H int {T n )). 

J n=0 ^ J J 
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The statistical averaging is carried out in (3) and (4) with 
respect to the zero-order density matrix of the conduction 
and impurity electrons. 
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The thermodynamic perturbation theory for Hi nt re- 
quires the generalization adequate for calculation of the 
statistical averages of the T - products of localized / - 
electron operators. This necessity appears for the rea- 
son that cannot be diagonalized with free / - electron 
operators. This Hamiltonian can be diagonalized by 
using the algebra of Hubbard' 28 ^ 30 ! transfer operators 
X mn — \m >< n| when the \m> with m = —1,0,1,2 
enumerates four states of the impurity atom: |0 > - is 
the empty or vacuum state with energy Eq = 0, the |1 > 
and |— 1 > or | |> and |!> are the states with one parti- 
cle with energy E a = e/ and spin a — ±1 and the state 
1 2 >=|f !> contains two / - electrons with opposite spins 



and the energy E 2 = U + 2e/. By using the relation 

.U = X° a +<?X° 2 , (7) 

we obtain the diagonalized form of the impurity Hamil- 
tonian 

2 2 

Hi = ]T E nX nn , Yl = L ( 8 ) 

n— — 1 n— — 1 

In zero order approximation, when we neglect the pro- 
cess of hybridization of the conduction and impurity elec- 
trons, the corresponding Green's functions have the form 
(u = u n = {2n + l)n/P) 



G° ff ,(k,k' \iuj) = <W<5kk 
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iuo — e(k) ' 
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where (a = —a) 

\ a (iLu) = ILU + Eq - E a 
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A^(jw) = + E w - E 2 , 

e -0^ + e -f3E 2 
n-a = 5? , 
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where n stands for (cr„, r„). The generalization for more 
complicate averages of type 

g° n (l,...,n | n + l,...,2n) = (-1)™ (Tf 1 ...f n f n+1 ..J 2n ) 
is straightforward, namely the right - hand part of this 
quantity will contain n! term of ordinary Wick type 
(chain diagrams) and also the different products of ir- 
reducible functions with the same total number of oper- 
ators. The full irreducible functions in 
g„(l, n\n + l,...,2n) also appears. For example 
5§(123|456) contains the contribution of 3! = 6 terms 
of ordinary Wick kind, then appear 9 terms of the form 

5°(l|4)^ 0)lr (23|56) and the last term is ^ 0)ir (123|456) 
The total number of terms is 16. In the case of 
(74(1 234 1 5678) there are 4! = 24 terms of ordinary Wick 

kind, the 72 terms of the type .g (l|5) ff °(2|6)^ 0)ir (34|78), 

then 18 terms of type 5 ^ 0)ir (12|56)^ 0)l, '(34|78), then 16 



In the case of / - electrons we formulate the identity 
which is just our GWT in this simple case: 
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terms of the form .g°(l|5).g^ 0)Jr (234|678) and finally one 
form gf )ir (1234|5678). The total number of terms is 131. 
The signs of all these contributions can be easily deter- 
mined. Thus the definition of the irreducible Green's 
functions or Kubo cumulants is just our GWT. In the ab- 
sence of Coulomb repulsion U all these irreducible func- 
tions arc equal to zero. When U 7^ they contain all 
the spin, charge and pairing fluctuations produced by 
the strong correlations. These definitions are the sim- 
plification of ones for Hubbard and other lattice models. 
The calculation of the simplest irreducible functions for 
example g^ lr (12\34) is rather cumbersome but straight- 
forward. It is necessary to find the values of chronological 
averages for 4! = 24 different orders of n , r 2 , t 3 and T4 
times and then to determine its Fourier representation 



(7/1/2/3/4)0 = <7/i/ 4 ) (7/2/3)0 - (7/i/ 3 ) (7/2/4)0 + (7/1/2/3/4)" - 

ff2 °(l,2|3,4) = .g (l|4).g (2|3)- ff (l|3) ff (2|4) + ff f ir (l,2|3,4), 
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^ 0),r [c7i, n; CT 2 , t 2 |ct 3 , t 3 ; cr 4 , t 4 ] = ^2 cxp(-zw 1 T 1 - iw 2 T 2 + iu 3 T 3 + i0J 4 T 4 ) X 

(12) 
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x 9 2 )lr [01 , iwi ; ct 2 , iw 2 1 03 , iuj 3 ; 04 , iw 4 ] ■ 
The Fourier representation conserves the frequencies 

g^ lr [ai,iui]a 2 ,iu 2 \(j 3 ,iuj 3 - 1 (j^iu)4\ = (35(u>i + u 2 - uj 3 - w 4 ) x 

(13) 

x g^ lr [ai, iuji, a 2 , iuj 2 \(7 3 , iuj 3 ; ct 4 , iui + iu 2 - iu 3 ]. 

I 

There is also the spin conservation u\ + a 2 = a 3 + 04. III. THERMODYNAMIC POTENTIAL 

Thus we have the rules to deal with chronological aver- 
ages of thermodynamic perturbation theory. First of all wc can determine the thermodynamic po- 
tential F of the system 

F = F -iln(C/(/3)) , 
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(14) 

F = -IlnZ -| Vln[l + exp(-/3e(k))], 

where Zo refers to the free impurity atom. The diagrams 
which determine the thermodynamic potential have not 
the external lines and are named vacuum. 

In Fig. I are shown the simplest vacuum connected di- 
agrams of the normal state. In the diagrams we shall 
depict the process of hybridization of C and / electrons. 
The zero order propagators of conduction and impu- 
rity electrons are represented by their solid and dashed 
lines correspondingly. These lines connect the crosses 
which depict the impurity states. To crosses are attached 
two arrows, one of which is ingoing and other outgoing. 
They depict the annihilation and creation electrons cor- 
respondingly. The index n means (ct„,t„) for impurity 
and (k n ,<7 n ,r„) for conduction electrons. The rectan- 
gles with 2n crosses depict the irreducible gn^ tr Green's 
functions. 

Besides the vacuum diagrams of fourth order shown on 
the Fig.l b) and c) there is also one disconnected dia- 
gram composed from two diagrams of the type Fig.l a) 
and containing additional factor 1/2!. Such situation is 
repeated in high order of perturbation theory and permit 
us to formulate linked cluster theorem. It has the form 

(f/(/3)) = exp(C/(/3)}^ (15) 

where (U{l3)) c Q contains only connected diagrams and is 
equal to zero when hybridization is absent. If we ad- 
mit the existence of the pairing mechanism of conduction 
electrons, thanks the hybridization, the paring mecha- 
nism appear also for impurity electrons. This mecha- 
nism results in appearance of the anomalous propagators 
of both kind of electrons. 

Fig. 2 shows some of the simplest connected anoma- 
lous vacuum diagrams. The anomalous propagators are 



depicted by the thin (solid and dashed) lines with two 
opposite directions at the end of them. 



IV. RENORMALIZED PROPAGATORS 

Now we shall consider the diagrammatical analysis of 
the perturbation series for renormalized propagators (3) 
and (4). The simplest contributions to such series are 
represented on the Figures 3 — 6. All such diagrams 
contain two external points with attached arrows deter- 
mined by the arguments of Green's functions and their 
kind. On the inner points of diagrams is supposed sum- 
mation on cr„,k„, and integration on r„. 

In the same second order approximation of perturba- 
tion theory the diagrams for impurity electron propaga- 
tors contain new diagrammatical elements namely the ir- 
reducible two particle Green's functions. These functions 
also can be normal or anomalous. The process of their 
renormalization will be not considered by us, supposing 
the necessity of renormalization only for the propagators. 

In Fig.5 the diagrams for impurity electron normal 
propagator are shown. 

The last two irreducible Green's functions of Fig.5 are 
anomalous ones because they contain non equal number 
of annihilation and creation /-operators enumerated in 
the left and right parts about the vertical bare corre- 
spondingly. Thanks the summation of infinite series dia- 
grams the renormalized normal and anomalous propaga- 
tors appear and now it is necessary to put equal to zero 
the source of electron pairs and simultaneously the bare 

f° and / together with anomalous irreducible Green's 
functions. The corresponding contribution to the anoma- 
lous impurity electron function /ercr'( r — T ') i s depicted 
on the Fig.6 

The final equations for renormalized functions it is 
more convenient to write down in Fourier representation 



G(k,cr,r|k',(T / ,r') = — ^ G aa > (k, k'| iuj) cxp -iu>(r - r ) 



F(k,<7,T|-kWy) - -^>^(k,-k'M 







exp 



-iuj(r — t ) 



The complete equations for the conduction electrons propagators have the form: 



G ffff '(k,k'M = <W<WG°(kM 



N 



^(G° CT (kMwMG°,(k'M 



- G°(kMjWMi^ v ,(-k'M) - F%(Miu:)^(-iu:)Fl,A-M^) - 

F CT -(kM/- CT ,MG°,(k'M), 



i^(k,-k'M = F^(k| iw)<w<W + 



N 



(Gl(k\iw)g aa ,(iw)F°,-,(k'\iw) + 



+ G°(kMiW,MG°,(-k'| - icj) + ^(kM. 9 ^Hc)G°,(-k'| - icj) - 

^kM7^,MF CT V<(k'M)- 



(16) 



(17) 



(a) 



(ft) 



FIG. 1: The simplest connected vacuum diagrams in normal state. The diagram a) is of second and b), c) of fourth order of 
the theory. Here V„ = V n /\/~N. 
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FIG. 2: The simplest vacuum anomalous diagrams. The diagrams a) and 6) are of second and c), d) and e) of fourth order of 
perturbation theory. 



These renormalized propagators are expressed through 
the full propagators g, f and / of impurity electrons. 
Now it is necessary to obtain the corresponding equa- 
tions for the full impurity electron propagators. Because 
the subsystem of /-electrons is strongly correlated we 
have to introduce the correlation functions Z aa i , Y a y and 
Yaa' which are represented by strong connected diagrams 
with irreducible Green's functions' 31-35 '. The process of 
renormalization of /-electron propagators is shown on the 
Figures 7 and 8, where the double dashed lines depict 
the full /-electron functions and the rectangles represent 

I 



the correlation functions K aa > = g® aa , + Z aa i , Y aa i and 
Y a a '■ The second equation we shall depict for anomalous 
propagator / of the impurity electrons (see Fig. 8). In 
both these equations the bare conduction electron prop- 
agators G®(k\iuj), F®-(k\iu) and 

F- a (— \n\ioj) play the role of mass operators for the /- 
electron propagators. It is easy to see that these func- 
tions participate in above equations being averaged on 
the Brillouin cell with matrix elements of hybridization. 
Therefore we define the new quantities 



^^^ 2 ^ kl G (k 1)CTl ,n|k 2 ,a 2! r 2 ) = l£ \Vu 1 \ 2 Gl lt72 (k\r 1 -r 2 ) = S ai(T2 G cri (r 1 -T2), 

kik2 ki 

i E V£T4V^(ki^i,n|k 2 ,a 2 ,r 2 ) = I^J^lUC-kiln-T,) = W 2 ^ 1CT1 (n -r 2 ), (18) 
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iVG< 2 >(k >CT ,r|k>V) = _ ^ V } V 2 G° 
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FIG. 3: Second order perturbation theory contribution for conduction electron normal propagator. 



NF {2 \k,a,r\ -k'.-a'.r') = — «-X-h — 



i G° V{ f V 2 * G° + F° VI g° V 2 * G ° 

l — « — X -« ► -X — ► ► — 1 1 ► — X ► — — ► x ► 



1 



F° Vi 7° V 2 F° 

-X — » — ■* X 



k, a, r 1 2 -k',-a',T 

FIG. 4: Second order perturbation theory contribution for conduction electron anomalous propagator. 

^ £ V^F ^,^^^,^) = l^l'^Ckiln -r 2 ) = ^^^-(n -r 2 ). 

kik 2 ki 

These definitions gives us the possibility to simplify the structure of equations for the /-electron propagators. By 
using these average bare propagators G° , F®— and and Fourier representation for r- variables we obtain 

A^iw) - Gi(-icj)[A CT (w)A ? (-icj) + y ffff (iw)Fy ff (iw)] 





d a (iuj) 


3 




b Fl a {iuj){K G {iio)K^{-iuj) - 






dcr(iuj) 




{Y aW (iuj) 


+ F°-(iLo)[(A a (iu)A w (-iLo) 


+ y £TF (iw)F FtT (iw)]} 


d a (iuj) 



d a {iu>) = {I - k a {iuj)G a a {iuj)){l - K^{-iLo)Gl{-iio)) + Y^{iw)F%{iuj) + 

+ Y aW (iuj)F% a (iu;) + F% a {iw)F%{i^)[Y^{w)Y^{w) + ^(iwJA^-iw)] + 

+ Gg:(-iw)Gj(iw)y o . ff (iw)y ffff (iw). (22) 

In the previous part of the paper we supposed the existence of pairing potential of conduction electrons with order 
parameter and with the bare propagators: 

iu + e(k) _ _o 



G ° (k| ^ } = M 2_^( k) ; ^( k M = = M 2_fl2 (k) ; ^( k ) = V^ 2 (k) + A2. (23) 

Now we shall discuss the case when the pairing potential is absent and the superconducting state appears simultane- 
ously with both subsystems as a consequence of the broken symmetry and phase transition. In this more simple case 
the renormalized conduction electron propagators have the form 

G^(k,k|tw) = Woa'tv( k l*w) + -^^( k Mw(*w)GS/(k|»a;) ) (24) 
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FIG. 5: The second order perturbation contribution for the impurity electron normal propagator. 
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FIG. 6: Anomalous impurity electron Green's function in the second order perturbation theory. 



G°(kM = (^-e(k))- 1 . 
The renormalized propagators of impurity electron in this case are: 

- GU-iw)[A a (iuj)A 7 (-iw) + Y a ^(iuj)Y^(iuj)} 



d a (iuj) 



d a (iuj) ' , " Jcrv 7 d a {iio) ' 
d a {iw) = (l-A .(iw)G^(iw))(l-A ff <-iw)G!2 f (-iw)) + 
+ G^(-w)G°(io;)y^(ia;)F^(ia;). 

I 



(25) 
(26) 

(27) 
(28) 

(29) 



The equation (24) has been established many years ago in the paper of Anderson^ by using the equation 
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FIG. 7: Dyson type equation for the normal propagator of impurity electrons. Double dashed lines depict full electron 
propagators. The arrows on them distinguish the normal and anomalous ones. The squares with attached arrows depict the 
correlated functions. On double repeated indices 1 and 2 is supposed summation by <r n and k n and integration by r„ . 
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FIG. 8: Dyson type equation for one of anomalous Green's functions of / -electrons. 



of motion of conduction electron operators. In this 
equation the propagator g a (iu)) has the role of t- matrix 



for non-spin-flip scattering. 
G CT (k, k'M 



By setting k = k' 



G ff (k,k'M = 



1 



\V k \ 2 g a (iu) 



(30) 



,~ ~ e(k) N(iu - e(k)) 2 

and considering the Lehmann spectral representation it 
is possible to conclude that the discontinuity of g a {E) 
across the real axis is pure imaginary! 8 ] 



9(T (E + iS) = \g a (E-i6)Y 



(31) 



The Green's function g a (iuj) has been known till now in 
approximate form as a result of special decoupling mech- 
anism used for equation of motion of quantum Green's 
functions. As is known in such decoupling approxima- 
tion some combinations of operators is taken off the av- 
erage value of product of operators and are replaced by 



their average values. After that truncation the Green's 
functions of low order remain. This approximation has 
been proposed by Bogoliubov, Tiablikov, Zubarev and 
Tserkovnikovl 21 " 24 ! and used by other authors' 2 " 14,18 !. 
The hybridization of conduction and impurity electrons 
causes the appearance of mixed Green's functions: 



G m (k,a,T\a ,t) = 
F m {k,a,T\a',T') = 
F m (-k, a, r\a', r') = 

and also 

G m (a,T\k,a',r') = 
F m (<j,r\-k,a\T') -- 

F">,T|k ;( T',T') = 

Let G m ^>(k\iu>), F ma - 



-(TC^(r)fAr'MP)) C , 
-{TC^{T)h,{r')Uml, (32) 

-- -(T.f a {r)C^{r')U{l3))l, 
-- -{TU(t)C^(t')U(0)}1$3) 
- -(TJ w (t)C^(t')U(P)) c . 

■(k\iuj) and Fmav' (k\iui) be the 
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FIG. 9: Schematic representation of the main approximations for the correlated functions The solid double lines with arrows 
depict the renormalized one-particle Green's functions of conduction electrons. The rectangles depict the irreducible Green's 
functions of impurity electrons. 










G° 






(0)ir 


— «-x— 




— x— «- 


(0)ir 
i/2 








G° 






X— «- 




— «-x— 




— x— *- 





(6) 



FIG. 10: Examples of the simplest ladder diagram for (/-function. 



Fourier representation of the first group of Green's func- fn the presence of superconducting pairing of conduction 
tions and G™ a , (k\iuj) , F™L, (— \l\iuj) and F Wr7 , (k\iuj) of the electrons we obtain the following results: 
second group. 
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G„w(kM = [Gl(k\iw)g^(iu)-F%(k\ioj)f^,(ioj)] 



F m ^,(k\iw) = ^[G°(k|^)/^M + ^V(kN^.9^Hc)], 

Y* r _ 

F m7f(T /(-k|iw) = -^\G^(-k\-iu;)f Wa ,(iuj) + F- a (-k\iu;)g aa ,(iuj) 



(34) 



For the second group of mixed propagators we have: 



G£,(kM = 



TV 



3 CTCT /(iw)G^(k|iw) - / (T ^(^)F-, (T ,(-k|^) 



0(-kM - -^[ w HF° ? ,(km + / a7 ( lw )GM-k|^)], 

F™,(k|i w ) = ^[/^MG°,(k|io;) + ^(-^ v ,(-k|^) 



(35) 



Now we multiply the system of operators (33) by V^/VN and sum after k, use the definitions (18) and suppose 
the paramagnetic phase of the system. Then we obtain: 

G"M = -1 J V k *(C(k|i W ) = Sff MGJH - /^(iwJ^H, 

F™ w (iw) = ^^VCF^(-Miw)=g a (iLj)F^(iw) + fM^)GU-^), (36) 
F™M = ^^y k *F™(k|^)=.g^(^)F° CT M+7^MG°M. 



When the superconducting state is established in the 
both subsystems simultaneously and the bare anomalous 
Green's functions of conduction electrons are equal to 
zero the above equations become more simple: 

G™M = g a (iuj)G a (iuj), 

F£{iu>) = /WMG&Hw), (37) 

For the second group of mixed functions in the same con- 
ditions we obtain: 

G mtT (iu>) = G a (iu)g a (iuj), 
F maW (icj) = G£(iw)/W(iw), (38) 
F mWa (icj) = G%{-iw)]^{ioj). 

V. APPROXIMATIONS 

In previous part of the paper we have formulated the 
Dyson type equations for the propagators of the system 
in general case of superconducting phase. These equa- 
tions contain the correlation functions which take into 
account charge, spin and pairing fluctuations and are de- 
termined by infinite sums of strong connected diagrams 
composed from irreducible Green's functions. The Dyson 



type equations for these correlated functions Z, Y and 
Y don't exist. Therefore to close the system of equa- 
tions and to determine the order parameters of the sys- 
tem state it is necessary to make some approximations. 
Our main approximations are determined by the dia- 
grams shown on the Fig. 9. 

Our approximations correspond to the summation of 
ladder diagrams in vertical direction shown on the Fig.10 
a). . We neglect the summation of ladder diagrams in the 
horizontal direction (see Fig.10 b)) 



VI. CONCLUSIONS 

The diagrammatic theory has been developed for one- 
site Anderson model in which strong correlations of im- 
purity electrons and their hybridization with conduction 
electrons is taken into account. 

The definition of irreducible Green's functions or 
Kubo cumulants is used as a generalized Wick theorem 
for strongly correlated subsystem of localized electrons. 
These irreducible functions contain all spin, charge and 
pairing fluctuations. On this base the linked cluster the- 
orem has been proved to determine the thermodynamic 
potential of the system and Dyson type equations were 
established for one-particle propagators of the electrons 
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of both subsystems. The main elements of these equations 
are the correlation functions Z aa i , Y a s' and Y-g a i which 
are composed from strong connected diagrams containing 
these irreducible Green's functions. 

The normal and superconducting phases are consid- 
ered. In the last case we examine the case when only the 
conduction electron subsystem has a pairing mechanism 
of superconductivity and when the superconductibility is 
established simultaneously in all the system as a result 
of broken symmetry. 
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